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A THEORETICAL STUDY OF THE LIFTING EFFICIENCY AT 
SUPERSONIC SPEEDS OF WINGS UTILIZING INDIRECT 
LIFT INDUCED BY VERTICAL SURFACES 
By Vernon J. Rossow 

SUMMARY 

A theoretical study is made of the possibility of improving the 
lifting efficiency of aircraft flying at supersonic speeds by using the 
indirect lift which can be induced on the wing by the fuselage shape 
and/or wing end plates. This investigation is carried out by analyzing 
a number of special models which have surfaces so located and inclined 
to the free stream that pressure shadows are cast over the wing surface. 
The major portion of this study is carried out using linearized or first- 
order theory. The method of calculating the pressures on each model is 
described and graphs of the final results are presented. 

The results indicate that in certain cases the lifting efficiency of 
a planar wing may be improved by favorable interference induced by lateral 
forces on the end plates or on the sides of the fuselage. At low values 
of reduced aspect ratio several of the airfoil systems studied are more 
efficient than is theoretically possible with planar wings. A configura- 
tion of comparable efficiency with a very high narrow fuselage is found 
for the higher values of reduced aspect ratio. In addition to improving 
the lifting efficiency of the aircraft, these indirect lift devices may 
also be used as control or stabilizing surfaces . 


INTRODUCTION 


Bodies or wings traveling at supersonic speeds cause disturbances that 
may extend for great distances without marked attenuation. Furthermore, a 
body traveling in the wave disturbance originated by another body may 
receive indirectly an aerodynamic force greater than that appearing on the 
body which originates the disturbance. ThiB property of supersonic dis- 
turbances leads to consideration of the possibility of' generating lift on 
a wing indirectly, or by interference. Effects of this type may possibly 
be used to increase the efficiency of the airplane by reducing the drag due 
to thickness or to lift. Either or both may be accomplished by favorable 
interference between the fuselage and the wing and in some cases between 
wing end plates and the wing. 
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The two-dimensional supersonic wind-tunnel nozzle is a classic example 
in which forces are induced indirectly. The pressure waves generated by 
the nozzle contour are cast across the flat wind-tunnel side walls which 
are alined with the stream direction. Since the side walls do not induce 
pressures, the side forces arise entirely from the indirect force induced 
by the nozzle contour. If a particular portion of this flow field is 
selected and placed in a supersonic stream, it could represent a wing 
(side-wall section) at zero angle of attack with lift induced on it by end 
plates (sections of the nozzle) which are deflected. It was the consider- 
ation of this analogy which led the author to the investigation reported 
herein . 

The use of end plates on wings at subsonic and supersonic speeds has 
been studied by numerous authors. For example, see references 1 through 4 , 
Reid, Hemke, Scholkemeir, and Gorgui., respectively. In general, these 
investigations considered the end plates as being fixed. Several of the 
studies included configurations wherein a portion of the end plate was 
movable. The movable parts were deflected in the same direction above and 
below the wing. The extension to "canted adjustable end plates" at sub- 
sonic speeds was investigated by Clements in reference 5. 1 The wing-end- 
plate combinations studied consisted of a straight wing with end plates. 

The parts of the end plates above and below the wing could be deflected 
independently to produce desired effects on the take-off and landing 
characteristics . 

The possibility of utilizing vertical surfaces on the fuselage to 
produce indirect lift at supersonic speeds has been studied by Ferri, 
Clarke, and Cassaccio in reference 6. The special configurations which 
they considered consist of a body with several flat sections that receive 
pressures from, and Induce pressures on, a wing. Two wing plan forms were 
considered, a delta wing and a particular tapered wing with the leading 
edges swept along the Mach lines . The trailing edge of the tapered wing 
was the hyperbolic intersection of the Mach cone from the body edge with 
the wing; The pressure fields were found by linearized theory and were 
integrated to obtain lift and drag. 

This report extends the study of interfering pressure fields at 
supersonic speeds by considering models that have swept wings, end plates, 
and specially shaped bodies . All these factors are varied to try to find 
an efficient lifting system. In a study of this type, the shape of the 
various components of the airfoil systems may be varied throughout wide 
limits so that the possible number of designs is unlimited. However, by 
choosing several models which are relatively easy to analyze. It is at 
least possible to discover which characteristics of the various models 
are desirable and to make improvements accordingly. 

1 The use of wing end plates deflected differentially for lift aug- 
mentation and roll control was proposed in a patent application by 
John A. Axelson of Ames Aeronautical Laboratory. 
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SYMBOLS 


A 

b 


°Di 

C P 

D 

G 

h 

i 

L 

m 

M 

P 

q 

Re 

S 

T 

u,v,v 

U 

x,y,z 


aspect ratio of wing, . ( . 8 P ari ) 

area 

one-half wing span 

chord length of wing 

lift coefficient 

coefficient of drag due to lift 

P - P~ 


pressure coefficient, 

' 4 00 

drag 

See equation ( 17 ) • 

height of rectangular part of fin 
(See sketch (j)-) 

lift 

ratio of tangent of leading- or trailing-edge angle to tangent 
of Mach angle (fj tan A) 

Mach number 

static pressure 

dynamic pressure 

real part of expression 

area of wing 

thrust on the fuselage 

velocity perturbations in x,y,z directions, respectively 

free-stream velocity 

rectangular coordinate system 

(The free stream is in the x direction.) 

See equation (l6) . 

See equation (6c) . 
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See equation (15) • 

See equation (23) and sketch (k) . 
angle of . attack 
^M 2 - 1 

deflection angle of fin or half-angle of cone 

optimum ratio of angle of attack of horizontal to vertical 
surface , ^ 


conical flow parameter, 
conical flow parameter, 


y + iz 


x + ‘J x 2 - y 2 
2e 


1 + e 2 

sweep angle of leading or trailing edge 


Subscripts 


free stream 

particular regions on wing 
body 

half-cone contribution - — 

cone interference on wing 
horizontal surface or wing 

quantities arising from interference of horizontal surface on 
body . 

quantities arising from interference of horizontal surface on 
cone 

quantities arising from interference of horizontal surface on 
itself 

quantities arising from interference of horizontal surface on 
vertical, surface 


leading-edge thrust 
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v-b quantities arising from, interference of vertical surface on body 

v-h quantities arising from interference of vertical surface on 

horizontal surface 

v-v quantities arising from vertical surface interfering on itself 

I, II, 

111 , 1 V, 

V 

Superscripts 


conical-flow solutions 

(See sketches (b) and (d).) 


( ) * first derivative with respect to deflection angle 

( )" second derivative with respect to deflection angle 


SCOPE 


The models that are studied may be divided into two general groups. 

The first group consists of combinations of flat vertical and horizontal 
surfaces of zero thickness without a fuselage. The second group includes 
a fuselage. A sketch of each of the models analyzed in this report is 
shown in figure 1 . 

The pressures acting on the various surfaces of the models are calcu- 
lated from conical-flow solutions derived by linear theory or by a super- 
position of several such solutions. A discussion of conical flow or of 
the method of finding the various solutions used, herein is not given 
because the subject is adequately covered in references 7 and 8 . The 
conical-flow solutions are taken directly from table A, 13a of reference 8, 2 
and applied to the problem at hand. The deflections of the various sur- 
faces are assumed small enough so that linear theory and the lin earized 
pressure relation. 



( 1 ) 


will apply with reasonable accuracy. Equation (l) is used throughout the 
analysis unless otherwise stated. In this way the lift and the drag coef- 
ficients and the drag factor of the various models are cal- 

culated for a range of reduced aspect ratios, gA. The coefficients are 

^The availability of this table avoids the necessity of finding each 
solution by linear theory and shortens the analysis considerably. 
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computed using both the plan- form area and the plan-form plus fin area, 
as reference area, thereby permitting a comparison to be made with stand- 
ard wings under conditions wherein the end plates or fins may or may not 
serve as necessary vertical-tail surfaces. When the vertical surface is 
considered as part of the fuselage, it is never included in the reference 
area. 


In order to simplify the calculations and the equations, the results 
are calculated for M = n/" 2. By use of the Prandtl-Glauert transformation, 
these results are extended to other Mach numbers, and shown in the various 
figures in a form that contains the Mach number as a parameter. 

The drag factor Cj^/^Cl 2 is used as the basis of comparison of the 

lifting efficiency of the various models. Its importance arises from its 
relation to the maximum lift-drag ratio expressed as 

fli\ = — — - ( 2 ) 

Wmax 2 n/c Do (C d ./c l 2) 

The coefficient Cp o expresses the drag of the whole airplane at zero 
lift, that is, when all the surfaces are undeflected. It includes the 
viscous and the thickness drag of the wing, fuselage, and auxiliary sur- 
faces. Effects of flow separation or possible effects of boundary- layer 
thickening on the various surfaces or surface junctures are not considered 
The coefficient of the pressure drag due to lift is and the lift 

coefficient is Cl* It Is assumed that Cl is proportional to the angle 
of attack and Cj^ to the square of the angle of attack, that is, a 

parabolic lift-drag polar. All of the models, except number 10, theo- 
retically have parabolic lift-drag polars . The advantage of using the 
drag factor Is that it is independent of angle of attack. It also affords 
a comparison of the lifting efficiency of the various airfoil combinations 
without fixing the thickness of the surfaces, ; the fuselage shape, or the 
viscous drag of the aircraft. 

It is necessary to establish reference models whereby the new lifting 
systems can be judged as worthy or unworthy of further consideration. The 
reference models chosen for the drag factor are the triangular wing, the 
constant -chord swept wing, and ..the theoretical minimum for slender wings 
of zero thickness. The results for these reference planar wings are shown 
on the various graphs as the daBhed lines . The parameters for the tri- 
angular wing are from reference 9> those for the constant-chord swept wing 
from reference 10, and those for the theoretical minimum from reference 11 

The reciprocity relations discussed in references 12 and 13 may be 
applied to all the models treated herein provided that in the reversed 
flow field the wing Is warped to carry the same lift. The drag and lift 
of the reversed flow model are the same as the model in forward motion 
and therefore the lifting efficiency will also be the same. The results 
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of this paper could then he extended to an equal number of models with 
modified wings which are flying in the opposite direction. This exten- 
sion will not "be made. 

The various fin-wing combinations and the method of calculating the 
pressures on each of them will now be described. The reader not inter- 
ested in the analysis of each model may go directly to the discussion of 
the results and the figures at the end of the report. The results of the 
calculations are presented in figures 2 through 7 • The drag factor which 
is used in the comparison is shown in figures 8 through 11 as a function 
of reduced aspect ratio. 

ANALYSIS OF AIRFOIL COMBINATIONS WITHOUT FUSELAGES 


In this section, several combinations of vertical and horizontal 
surfaces are described and analyzed. The general procedure followed in 
each case is: 

1. Model is described. 

2. Wing and fin flow fields are divided into the various conical 

flows which are superimposed to determine the pressure on 
each element of surface. 

3 . Expression for pressure in each conical flow field is found from 

the general formulas of reference 8 . 

4. Expressions are combined and integrated for lift and drag 

coefficients . 


Model 1 - Triangular Fins With Unswept Leading 
Edges and Sonic Trailing Edges 


Description .- The first airfoil 
combination to be studied consists of 
a wing swept along the Mach lines 
with end plates or fins as shown in 
sketch (a) . The fins have unswept ^ 

leading edges and the trailing edges ; — 
are swept along the Mach line. The 
wing is parallel to the free stream 
and therefore does not generate a pres- 
sure field directly. The portion of 
the fin extending above the horizontal 
surface is deflected so that an expan- 
sion field or negative pressure region 
is produced on the top of the horizon- 
tal surface or wing. The portion of 



Sketch (a) 
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the fin extending below the wing is deflected so that a compression field 
or positive pressure region is produced on the bottom of the wing. The 
wing serves as a barrier between these opposite pressure fields. The 
deflections of both the upper and the lower parts of each end plate are 
always taken as equal and are considered, positive when the lift generated 
by each is positive. 


Lift.- On the wing of this first model, "the pressure distribution is 
of the two-dimensional or Ackeret type ahead "of the Mach cone from the fin 
tip. Inside this Mach cone the pressures take on a more complex form. 


As pointed out previously, all edges of 'this configuration are either 
sonic or supersonic and therefore the flow field between the fins on top 
of the wing may be treated independently of the flow between the bottom 
fins. Also, since the upper and lower fins are deflected the same amount, 
the lift contribution induced by each may be “Considered equal within the 
linear approximation used in this analysis. The indirect lifting pressures 
within the Mach cone from the fin tip for the top side of the wing are cal- 
culated by superposition of several basic conical- flow solutions as illus- 
trated in sketch (b) . The values of u and w shown in the various regions 
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^ Plane 
of 

Symmetry 

i n m 

Sketch (b) 

are the boundary conditions in the plane of the surface, z - 0. Addition 
of these flow solutions yields the solution together with the appropriate 
boundary conditions for the upper half of the fin and its image. 




Sketch (c) 


The resulting fin shape shown on the 
right of the sketch has a plane of sym- 
metry which may be taken as a boundary. 

The wing is assumed to lie in this plane 
and serves as a reflecting plane. One 
half of the resultant flow field shown in 
sketch (b) is equivalent to the flow field 
between the fin and wing on the top (or 
bottom side," depending on which part of 
the flow field was chosen) as shown in 
sketch (c) . The velocities In the .flow 
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field inside the Mach cone from the fin tip which are induced by the 
conical surfaces shown in sketch (b) are found from table A, 13a of 
reference 8 . The axial velocity for the first two parts are given by 
number 9 of this table as 


(3) 

mi — > i.o 

mg oo 


u-r = Re 


2m 2 US 


it "•/ 


tan 


-i /(m g - 1 )(t - m x ) 


mg 


2-1 


V (m 2 - mi) (l - t) 


where m x and m 2 are the relative sweep angles of the trailing and. 
leading edges , respectively , and 

m = (3 tan A (4a) 


T 


2 e 

1 + e 2 


(to) 


e 


y+ 12 _ 

X + n/x 2 - y2 _ 22 


(to) 


After substitution of equations (to) and (to) into equation ( 3 ) and finding 
the real part at the above limits, the perturbation velocity in the x 
direction is. 


U I = 


US 

it 


tan -1 z y - + 2 tan -1 

x2 - Z2 + X •>/ x2 - y2 - z2 


X - y 


+ n/"x^ 


- Z 2 

(5a) 


) 


This expression is valid within the Mach cone from the fin tip; that is, 
x 2 > y2 + Z 2. Ahead of the Mach cone, x 2 < y 2 + z 2 , the square root 
becomes imaginary and equation (5a) is not valid. However, the velocity 
ahead of the Mach cone is known, to be the same as for two-dimensional flow, 

Uj = US for z < x < -v/y 2 + z 2 (5b) 

On the surface. of the wing, y = c, inside of the Mach cone from the fin 
tip the pressure coefficient is 



since u-q is the same as uj in the y = c plane (i.e., the reflection 
plane or the plane of the wing) and "in/ 11 - 8 - 
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The pressures on the wing both inside of and forward of the Mach cone 
from the fin tip are now defined . The lift on the wing in these regions 
may he divided into two parts, A and B, as shown in sketch (c). Part A 
is the lift induced on the wing surface between z = 0 and z = h and 
part B is the lift on the wing between z = b and b + (c/2). Region B 
is noted to be beyond the center of the whole configuration. On this 
part, lift is induced by fins at both wing tips, resulting in a dcrubling- 
up effect over the center sections. These parts may be written separately 
as, - 



vc 2 + z2 ) dz + 46 -jp 

where 


(6b) 


5 = (2b + c) g - c 2 
2(2b + c ) 


(6c) 



f 


r 


The total lift coefficient for the entire wing Is the sum of equations (6a) and (6b), where the 
coefficient Is based on plan area plus fin area, 2bc + 2c 2 . The integration with respect to x 
can be carried out easily. So, equations (6a) and (6b) become 


*b 

% = cH) l {' tan-'Vf - « + (. + c)[tan-l 

. c + z + •/2cz j , c - z 1 . 

z Jn • — — — + c sin 1 ■■ ■ - dz 

v c s + z 2 c z + z 2 J 


zc „ , z I 

== + 2 tan pm + 

c2 + 2cz + (c + z) v 2cz z +V2czj 


r c 2 + z 2 


°% ’ 5^ (j + if { (2b 


+ c - z) tan -1 = — + 

L (2b + c - z) 2 - z 2 + (2b + c - z)*J (2b + c - z) 2 - c 2 - z 2 


2 tan" 


(2b - z) + *J(2b + c - z) 2 - c 2 - z 2 


-*>/c 2 + Z2+ z tan-^( ? b+ c : z )l,:, cS .r . 


z In 2 +tJ W + c " z ) , 2 - c2 - z ! + c Bln -X c(c + 2b - z - c) - z 2 + *c _ 

n/c 2 + z 2 -7 C 2 + z 2 (2b + c - z - c) 2 


} a j) 


n(2b + c - z - 4 c z + z 2 ) 


The Integration over z is carried out graphically. 

Since the wing has a constant chord, the aspect ratio is, 

A - Tf (8) 

As pointed out above, these results are for M = V2. However, by application of the Prandtl- 
Glauert transformation the results of equations (7) and (8) can be extended to other Mach numbers. 
This transformation, of course, changes the proportions of all the airfoils to fit the changed 
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geometry of the Mach waves. All the results of this paper will be extended 
in this manner without further reference to the Prandtl-Glauert rule. 


Drag .- The horizontal surface on which the lift is acting is at zero 
angle of attack and of zero thickness. Therefore, it does not contribute 
to the pressure drag. The only drag arising is that of the deflected fins. 
It is to be noted that the inner side of all the fins is not influenced by 
any of the other surfaces. The pressure distribution on these surfaces is 
therefore the two-dimensional type, Cp = ±26. 


On the outer sides of the fins, velocities are induced by both the 
upper and the lower fins over a portion of the fin surfaces. Since they 
are deflected differentially, the drag is considerably less in this region 
of mutual influence. The drag of the outside of the fins is computed by 
superposition of the several conical flows as shown in sketch (d). The 
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leading edge of solution V is alined with, the trailing edge of solution I 
so that the downwash from IV is completely canceled. Solution I has 
already teen discussed and the streamwise velocity perturbation is given 
ty equation (5a) • Solutions IV and V are given ty numbers 1 and 9 of 
table A, 13 a of reference 8. Since Uy does not influence the pressures 
on the fins, it will not be considered at this time. 


u IV = Re 


mUS 


it \/m 2 - 1 


cos 


-1 


1 - mT 


m - t 


m — > 00 


(9a) 


Upon allowing m to proceed to the limit and taking the real part 


2 ® (2 - tan - 1 Z JL - + 

\2 x 2 - z 2 + x s/ x 2 - y 2 - z 2 

2 tan -1 — ^ ( 9 b) 

x - z + x 2 - y 2 - z 2 / 


u iv| z=0 “ IT (| + sin_1 x) < 9c) 

Therefore, on the surface of the fin (z = 0) the velocity perturbations 
are. 


u = +U 8 ahead of mutual Influence region (10a) 
u = ± 222. sin ” 1 ^ inside of mutual Influence region (lOb) 

By use of equations (l) and (10) the integration of the pressure over 
the fin surfaces can be carried out to yield the drag coefficient based on 
plan area plus fin area. 


_ (Vtc + 2)5 g c g = 3.2732S g c a 
be ■+ c 2 be + c 2 


( 11 ) 


The drag factor, C^/ pC-^ 2 , is calculated from the results of 
equations ( 7 ) and (ll) . 
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Model 2 - Triangular Fins With Unswept Leading Edges and 
Sonic Trailing Edges; Wing With Curved Trailing Edge 



Description .- On the rear portion 
of the wing of Model 1, shown in 
sketch (a), there is a region of nega- 
tive lift. In this section, the wing 
plan form is modified as illustrated 
in sketch (e) so that the trailing edge 
lies along the line where the local lift 
becomes negative. The remainder of this 
airfoil system is the same as Model 1. 

Lift .- The equation defining the 
pressures On the horizontal surface of 
Model 2 is the same as for number 1, 
that is, equation (5) • From this equa- 
tion it is possible to find the expres- 
sion for the line on the wing along 
which Cp = 0. 


Sketch (e) 


Cp = 0 = ±2S T 


_4S 

rt 


tan 


-l 


zc 


X 2 - Z 2 + X \lx 2 - c 2 - z 2 


2 


tan -1 


x 


c + 



Through the use of trigonometric relations and after some algebraic manipu- 
lations, the equation for the line on the wing along which the pressure 
is zero is 


= nTx W 


xc 


or 


x = 


c + n/ c 2 + 4z 2 


(12a) 

(12b) 


The plan area of the wing is now found as 

*b, 

'o 


S = 2 J' ■7 c 2 + hz 2 - z^dz 


o , b f r? 1,-u'o , ^ 2b + C2 + lft} ^ 

S = cb + — v c 2 + 4b 2 - b 2 + — r> in 


( 13 ) 



I 


1 


w . 1 

a 


g" ( ( ' + ' 


IXCtylVtl U 


Region A 


The lift is found by integrating equation ( 5 ) over the wing 
surface. As before, the lift is divided into two parts, A and B, 
as indicated in sketch (f) . The coefficients are, after inte- 
gration with respect to x 


I— c— 1 


■Intersection of Mach 
Cone with Wing 


-x=-|-(c+j 6 i + 4 P) 


Sketch (f) 


, _ 16& 
irS 




■ hz 2 - c rtz . Itc . „ c + -VcS + hzs. +v2c(vc2 + kzz - c) , 

oZ " ~o~ + ~o~ + 2 in 7== : : u + 

2c d d 2 n/c 2 + z 2 


c sin-l c(^ C g + 4z 2 - c) - 2zs ] dz 
■s/c2 + Z 2 ( s/ c2 + 4 z2 - c) J 


(lM 


Ss’ff + Tjf 


+ 2 tan" 1 


*m E “ Z 2 + - c2 - z2 


3%! - C + V Xjn 2 „ c 2 - z 2 • 


jtXjn + ~ + z, tan' 


■1 / X * S ~ e2 - z - + z u x m + ^ 2 ~-^jg + 0 Bin -x - c) - zs 


' C2 + Z2 


c 2 + z 2 (xm - c) 


M 




where 


lfb(4b 2 - c 2 ) + 2*/2bcs/2b 2 + c 2 

— 


= f c + 


(16) 
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and G is the lift over region B when a two-dimensional pressure dis- 
tribution is assumed to exist. 


G = 25 


be 2 


4b - c 


- (whj + + + - tstzJ 


Wc 2 + 4b g + £f j n 4b c 
2 4 


- 8b g + (4b - c) 


y- ♦ <■> - 


(4b - c)(s/ c 2 + 4b 2 - 2b) 


(17) 


The aspect ratio for Model 2 is, 



where S is given by equation (13)* 


( 18 ) 


Drag . - The pressure drag of the fins of Model 2 is the same as for 
Model 1. The coefficient is different because of the change in the plan 
area. So 




v-v 


^ + 2^6 2 c 2 

S + 2c 2 


6. 5465 2 c 2 

S + 2c 2 


(19) 


Model 3 - Triangular Fins With Sonic Leading 
Edges and Unswept Trailing Edges 



Description .- The horizontal sur- 
face of this airfoil combination is 
the same as number 1. The fin shape 
differs in that the leading edge is 
swept along the Mach line and the 
trailing edge is unswept as shown in 
sketch (g). 

Lift.- Since all the edges of 
this configuration are either sonic or 
supersonic, the upper and lower sur- 
faces of the wing may once again be 
treated independently. 


JL 

IT 


Sketch (g) 
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The boundary condition on the 
fin and on the wing can he satisfied 
by the selection of a plane of sym- 
metry from a particular conical flow 
as illustrated in sketch (h) . The 
Mach cone from the fin tip does not 
intersect the horizontal surface. 
Therefore, the streanrwise velocity 
perturbation in the y = 0 plane 
yields the pressure on the wing. 

From number 3 of table A, 13a 
of reference 8, the streamwise veloc- 
ity perturbation for the conical fin 
shown in sketch (h) is 



= Re 


mU6 


it •J 


m-= 


( cos 

- 1 \ 


■1 1 - IT 


m - t 


+ COB -1 


1 + mT 
m + t 


m 


After m has proceeded to the limit and the real part is taken through 
the use of equation (4), 


2U6 (x 2 - z 2 + x 

"v " — 


x a - y g - z 2 )(x 2 - y 2 + x J x a - y 2 
(x 2 - y 2 + x ■/x 2 - y 2 - z 2 ) a + z2y 2 


- zg) - z2yg 


( 20 ) 


On the top of the wing, y = 0, 


1 = _ _ 4S si x 2 - z 2 

'P " TJ jt x 


( 21 ) 


As was done with the previous models, the 
lift is divided into two parts, A and B, 
part A being direct lift and B the overlap 
lift (sketch (i)). The lift for these two 
parts is integrated and found as 



H — c — H 


Sketch (i) 
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8S 

rt(bc + c 2 ) 


-fz be + e 2 




sin” 1 


_s_V 

b + c/ 


(2b + c) 2 cot- 1 yi^-e 


(22a) 


Cr — 


85 


«(bc + c 2 ) 



2bc + c 2 ) 


3 / 2 


itg- l- J 3 - s/gbc + c 2 - 
6 




Drag .- The pressure on the inside of the fin surfaces can be calcu- 
lated from equation (2l) by setting z = 0. On the outside surfaces, the 
upper and lower fins interfere on each other in such a way as to reduce 
their drag. The pressures on these surfaces could be found and integrated 
for the drag. This procedure is not necessary because the drag can be 
found directly from the results of the previous sections through the use 
of reverse-flow theorems. Inspection of sketches (a) and (g) shows that 
the fins are identical If the flow direction over one of them Is reversed. 
Within the linearized approximation being used in this report, the drag 
of the two fin plan forms Is the same in forward and reverse flow as 
demonstrated in references 12 and 13 . Therefore, the drag coefficient for 
this wing-fin combination Is also given by equation (ll). 


Model 4 - Trapezoidal Fins on a Swept Wing 



Sketch (j) 


Description .- The three previous 
models have fins which are only of the 
order of one chord length high. Con- 
sider now the case when the triangular 
fin is raised and a rectangular section 
Is added as shown in sketch (j). The 
upper and lower fins are deflected 
differentially to give a positive lift. 
All edges of the fins are unswept 
except the fin tip, which is swept 
along the Mach line. 


Lift and drag .- The fins are 
assumed to be high enough so that the 
Mach cones from the tips do not inter- 
sect the wing. The minimum height, h, 
CONFIDENCES** 
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of the rectangular section which affords this condition is found by seeking 
the point where the Mach cone from the fin tip first intersects the wing. 
The equation for the intersection of the Mach cone with the plane of the 
wing is. 


x 2 = (c + h) 2 + z 2 


when the Mach cone first makes contact 
with the wing trailing edge, x = c + z 
(sketch (k)). Substituting this in 
the above equation 

z - h + |£ = z h ( 23 ) 

If the fin is high enough so that 
z^ > b and h > c, the lift and drag coef- 
ficients based on plan form plus fin area 
are, 

48b c + 48 

4- 

^Lv-h be + c 2 + 2hc 



Sketch (k) 


(24a) 


4S 2 (2hc + c 2 ) - ~ B 2 c 2 
CDy-v be + c 2 + 2hc 


(24b) 


The second term in the numerator of the lift coefficient Is the lift con- 
tribution from the doubling-up effect at and near the center of the wing. 
The second term, in the numerator of the drag coefficient is the reduction 
in drag due to the mutual interference of pressure fields of the upper 
and lower fins in the region indicated by the shaded area In sketch ( j ) . 


Model 5 - Wing and Fins in Horizontal Plane 


Description .- The question may be asked 
whether it might not be profitable to lean 
the fins outward with respect to the wing as 
i llu strated in a front view of sketch (l) . 



Front View 


Sketch ( l ) 
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o 


Plan 


Front 



Side 


The general case is difficult to eval- 
uate because of the superpositions 
involved. However , the limiting case 
of the fins and wing in the horizontal 
plane, sketch (m) , is simple to ana- 
lyze. The fins (one on each wing tip) 
are deflected so that they not only 
induce lift on the wing, which is at 
a = 0, but carry lift themselves. 

Lift and drag .- The pressures on 
the fins and wing are given by equa- 
tion (9c). The integration of these 
pressures for lift and drag in the 
various regions is carried out. 


Sketch (m) 


C L - 


— {o + 

be + — 


C 2 ) + - 
’ Jt 


b 2 sin -1 — - — + 2(c 2 + 2bc + 2b 2 ) 
d + c 


(f - tan - 1 | (b 2 + be + c 2 ) - b ‘J 2b c + C 2 J- 


(25) 




1.6375 a c 2 

j c 2 
be + ~ 


( 26 ) 


The aspect ratio is given by 


A = (2b + 2c ) 2 
(2b + c) c 


(27) 


Optimum Deflection of Wing 


Description .- In the airfoil combinations studied so far, the wing 
has been undeflected. It is therefore not being utilized to the best 
advantage. In order to obtain the highest efficiency from the wing, it 
should be deflected so that the drag factor C r j 1 /j3C L 2 is brought to a 
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minimum - This minimum is derived "by finding an optimum ratio, 0 = a/8 , 
between the wing angle of attack and the fin deflection. When both the 
wing and fins are deflected, the drag factor may be written as 



The primes in equation ( 28 ) signify derivatives with respect to the deflec- 
tion angle of the various surfaces. Note that there is a contribution to 
the drag arising from lift induced by the fins on the now deflected hori- 
zontal surface. 3 Differentiating equation ( 28 ) with respect to a while 
holding all the coefficients constant and then setting the result equal 
to zero furnishes the maximums and minimums of the Cp^/ Cp 2 versuB a 

curve. Two values of a are found: When Cp^/Cp 2 is a maximum, that 
is, not an optimum. 


a 


<Vh 

°Lh-h ' 


% 

C L 2 


OO 


and when Cp^/Cp 2 


is a minimum, an optimum. 



The first value of 0 defines the ratio in which the fins and wing should 
be deflected to yield zero lift; that is, it defines how the wing is to 
lift an amount equal and opposite to the lift induced by the fins. There- 
fore, the model has drag without lift and the drag factor becomes infinite. 

The optimum value of 0 defines the ratio in which the fins and wing 
should be deflected to give the most lift for 'the drag. In general, 0 
is positive, indicating that both the fins and the wing lift contributions 
are in the same direction. Hereafter, the term 0 will be used to denote 
only the value which yields the optimum ratio. 

3 Considerations of this type must be taken into account when thick- 
ness of the various surfaces is considered. 
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Equation (29) could be applied to all of the airfoil comb in ations 
studied so far except number 5. It is applied to Models 1, 2, and 3, 
but not to 4. Before this is done, it is first necessary to find C 

the lift on the wing when it is deflected, and the drag C-n 

■^h-v 

on the fins by the wing. Since all the edges are either sonic or super- 
sonic , none of the surfaces have leading-edge suction and therefore C-q 
is zero for these three cases. ^ 


L h-h' 
induced 



Lift of finned wing .- The pres- 
sures on the wing of Models 1 arid 3 
are found from a superposition of 
several conical flow fields. It will 
be assumed here that the fins are 
undeflected and serve merely as bar- 
riers or reflection planes. It then 
remains to choose a symmetrical flow 
field from which the parts can be 
dissected. A flow field of this type 
is given by number 3 of table A, 13a 
of reference 8. The streamwise per- 
turbation velocity has already been 
given in equation (20) . The conical- 
flow plan form and the part of the 
flow field to be used is shown in 
sketch (n) . The manner in which the 
wing flow field is constructed is shown 
in sketch (o) . The streamwise velocity 
perturbation on the wing is found by 
setting z = 0 in equation (20) . 




By use of equations (l) and (30) , the lift of the finned wing is found by 
integrating the pressures over the various surfaces shown in sketch (o). 

c _ 8a, 2(b + c) s/2bc + c g - 2 c g 
^h-h 3*: be + c 2 ' 


(31a) 



for Models 1 and 3 • For loth models the drag is simply. 


C-n = <xCt 
"h-h 4i-h 


(31b) 


The plan, form of the wing for Model 2 is different from 1 and 3 hut 
the same expression for the pressure is used. The limits of integration 
for the lift in this case are not so simple as for Models 1 and 3. There- 
fore, the spanwise integration is carried out numerically. The same 
superposition technique shown in sketch (o) is used for the wing with the 
curved trailing edges (still supersonic), and the lift is found as 


JL h-h 


16a 

TtS 



c + 




4y £ 


- y 2 iy + 


4bz 


p 4b-c 1 

c + v c 2 + 4(2b - y) 2 

2 

1 J 

L 2 



- y 2 dy + 



1 


c +«Jc 2 -t 4(b - y) 2 


1 2 

-b - 
_ 



(32) 


where S is defined by equation ( 13) • The values computed by equa- 
tion ( 32 ) are used in the optimum case of airfoil system 2 and are not 
shown in figure 2. 


Interference of wing on fins .- Heretofore, when the lift from one 
surface was considered, the other surface was assumed to be* undeflected. 
If both are deflected simultaneously, drag contributions will arise due 
to one surface interfering on another. As already pointed out, the lift 


induced by the fins on the wing contributes to the drag as 
other cross inteference arises from the wing 
interfering on the deflected fins, desig- 
nated in equation (28) as Of! 

h-v 

The pressures induced on the fins by 
the wing are given by equation (2l) . The 
region on the inside of the fins over which 
the pressures act is shown for the various 
configurations in sketch (p) . After inte- '- x 
gration of equation (2l) over these areas, 
the interfering drag is Models 1 and 2 



5ffC W 


The 



Model 3 


Sketch (p) 
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Model 1: 


Model 2: 


20 - l)a8c= 


_ . 0 *5k6aSc z 

J h-v be + c2 be + c 2 


~ ^ a ® c2 ^ 


093aSc 2 


J h-v 


Model 3: 


"D 


h-v 


a&c 2 
be + e 2 


(33a) 


( 3313 ) 


(33c) 


With all the factors in equation (28) now known, the optimum relative 
deflection angles and the drag factor are found. 


ANALYSIS OF AIRFOIL COMBINATIONS WITH FUSELAGES 



The same principles that were 
applied to the airfoil systems in the 
last section may be used in contouring 
the fuselages of aircraft. Several of 
these models are now described and 
analyzed. The fuselage is treated as 
thin in cross section and, in general, 
semi-infinite in length. These thin 
bodies which are tall in comparison 
with their, width, are thought of as 
perturbations on the fuselage shape 
to be used in a design, as indicated 
in sketch (q). If the fuselage is 
slender (maximum diameter approximately 
equal to height of perturbation body) , 
this method of treating the problem 
is not believed to affect greatly the 
efficiency of the configuration. 


Sketch (q) . 
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Model 6 - Triangular Fins With Unswept Leading 
Edges and Sonic Trailing Edges 


Description .- An examination of the 
streamline paths over the wing of Model 1 as 
indicated in plan view of sketch (r) , reveals 
one criterion for contouring the fuselage. 

If a slender body is placed along the center 
line of the model, the part on top of the wing 
must increase and the part on the bottom must 
decrease in width to follow the streamlines. 
These lateral perturbations which would be 
required on a conventional body are illus- 
trated in sketch (s). The pressures induced 
by the fins exert a thrust on this new con- 
tour. The paths of the streamlines and thrust 
on the body will now be calculated. 



Thrust on fuselage due to fin 
deflection .- For simplicity of cal- 
culation the body is assumed to lie 
in the z = b plane. The local 
pressures on the fuselage are given 
by equations ( 5 ) . Since u-j- ^ u-j-j- 
at all points on the body, a separate 
expression must be written for the 
pressure in each of the three regions 
shown in sketch (t) . Region 1 
includes that part of the fuselage 
which feels the full induced pressure 
of the fins but does not realize that 
fins are finite in extent. In region 2, 
the pressures are decreasing in the stream 
direction because the tips of the fins are 
relieving the pressure more and more. The 
boundary of region 2 is defined by the inter 
section of the fin-tip Mach cones on the 
fuselage and the reflection of these Mach 
cones from the wing surface. In region 3 
the pressures decrease still further. This 
reflection effect is analyzed by placing an 
image fin below the wing. 



-Reflected 
Mach Cone 
, Intersection 


(Dv 


(2) 


(3) 


Original 
Mach Cone 
Intersection 

Sketch (t) 


The pressure induced by the fins in each of the three regions on the 
upper body are 
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Region 1: 


Region 2: 


C p = -26 


(34a) 


Cp - - -f f tan -1 


^L===m====+ 2 tan -1 

x 2 -b 2 +x \l x^-y^-b 2 x-y+ \l x 2 


-y 2 -b 2 ^) 


Region 3: 


P _ 28 
Cp J 


(34b) 


tan" 


by 


x 2 -b 2 +x V x 2 -y 2 -b 2 


,+2 tan 


— z 


tan 


-l 


b(2c-y) 


x 2 -b 2 +x */ x 2 -(2c-y) 2 -b 2 


+ 2 tan 


-i 


x-y+ n/ x 2 -y 2 -b 2 
b 


x+y-2c+ *</x 2 -(2c-y) 2 


-b 2 . 


+ 26 


(34c) 


The thrust on both sides of the upper fuselage may be written approxi- 
mately as 


! - */T 

J o '•'b 


c pb+c 


w 


Cp ^ dx dy 


(35) 


where w/u = dz/dx and the equations for the pressure and slope change 
as the integration proceeds over the various regions . 

The slope of the body, or streamlines in the z = b plane, is found 
by substituting the proper limits into solution number 9 of table A, 13a 
of reference 8. 


Vt = Re 2UB rg i ~ m i Binh -i / mi - t .^-i fc-iaJd+r) ' 

* L* miUg+l) -/(l+n+Jr v (m 2 +l) (i,-r) J 


m x l 

mg — ;> OO 


(36) 


After the limits are inserted and the real part is taken through the use 
of equation (4) 

(t) = i (| ln X+ ^ x 2 -y g - z2 . 4.^-1 x-y+z+^xa-ye-za ^^-1 x-y-z+ sTxg-yi^i N 
^ x 2 -y 2 -z 2 x+y+z+ v x 2 -y 2 - z 2 x+y— z+ \/ x 2 -v2— z2 / 
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The slope of the hody in each of the three regions shown in sketch (t) 
is then 


Region 1: 


i - 8 


( 38 a) 


Region 2: 


S'® 


(38b) 


z=h 


Region 3: 


w _ / w 

U " lu. 


Ar\ t 8 f 1 ^ x+ ^ x 2 -(2c-y) g -b 2 ] tan _i x+y-2c+b+ N/x g -(2c-y) 2 -b 2 

\V-r ^ j_2 x _ x 2 -(2c-y) 2 -b 2 x-y+2c+b+ x 2 -(2c-y) 2 -b 2 

iz=b 


tan' 


-1 x+y-2c-b+ s/x 2 -( 2c-y) S -b 2 
x-y+2c-b+ *7 x 2 -(2c-y) 2 -b 2 


( 38 c) 


The complexity of equations (34) and ( 38 ) necessitates the integra- 
tion of equation (35) by graphical means. The local pressure, slope, and 
hence thrust are calculated at various points on the fuselage. These 
local thrust values are then plotted and integrated for total thrust. 

The thrust on the lower part of the fuselage is the same as that on the 
upper part so only one integration is necessary. Typical values of slope 
pressure, and thrust for the upper body are shown in figure 5 f° r b/c = 
Note that negative values arise in region 3- All integrations are stopped 
when the pressure or slope on the body changes sign. The final integrated 
values of the thrust are shown in figure 6. 

The equation for the slope of the fuselage can be integrated to yield 
the coor din ates of the fuselage shape. Since w/u = dz/dx and because 
the z variations are assumed to be small, z may be set equal to b for 
the integration of equation (37)* 

Region 1, x < '/yiTb 2 : 


z = 8 ( x-b ) 


(39a) 


*\ H 
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Region 2, x > "7 y 2 +b i 
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(39b) 


The coordinates of the body in region 3 are given by an equation which is 
equivalent to equation (39b) plus similar terms to account for the 
reflected Mach cone interference. 

Thrust on fuselage due to fin and wing deflection .- When the wing 
is deflected in combination with the fins, pressures are induced by the 
wing which also exert a thrust on the fuselage. The thrust on the per- 
turbation body is expressed as the integral of the product of the slope 
induced by the fins and the pressure induced by the wing over the fuse- 
lage sides. This thrust augments that due to the fin interference. 

From equation (20) the pressure on the fuselage induced by the wing 
is 


°P " 


-Ifa, ,x g -(y-c) 2 +x , 7x g -(y-c) g -b g . _x 2 -b 2 +x : *J x 2 -(y-c) g -b a . - (y-c^b 2 


it 


Jx 2 -b 2 +x *7 x 2 -(y-c) g -b 2 2 + (y-c) g b 2 


(^ 


where the change in coordinates has been incorporated so that equa- 
tions ( 38 ) and (to) correspond to the same coordinate system. As before, 
the integration is carried out graphically. The results of this integra- 
tion are shown in. figure 6. The thrust coefficient based on plan form 
plus fin area is shown in figure ^(b) along with the drag coefficient 
of some of the other models discussed. The thrust coefficient for the 
fin interference on the body and for the wing interference on the body 
is represented by ^ and respectively . 

The lift and drag for Model 6 is given by equations (7a) and (ll) 
when only the fins are deflected. The overlap" lift given by equation (7b) 
has been eliminated. The optimum combination is found by equations ( 28 ) 
and ( 29 ) along with ( 31 ) for the horizontal-surface values. 
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Throughout the analysis of the thrust of this wing-body combination 
it is assumed that the slope of the sides of the perturbation fuselage 
deviate from the vertical by a negligible amount . 


Model 7 - Triangular Fins With Unswept Leading Edges and 
Sonic Trailing Edges; Wing With Curved Trailing Edge 


Description . - This model, 
sketch (u) , is the same as Model 6 
with the exception that that portion 
of the wing has been removed which 
contributes negative lift, as already 
discussed for Model 2. The forces on 
the various components are found from 
the relations already developed. The 
thrust and shape of this fuselage are 
the same as Model 6 because the inte- 
gration is stopped when the pressure 
becomes negative, or at the trailing 
edge in this case. The lift is found 
from the results of Model 2. 



Model 8 - Positive and Negative Wedge Superimposed on Fuselage 


Description .- The pressure waves 
were cast out over the wing by the 
fin in the last group. The same 
effect could be obtained by con- 
touring the fuselage to do this. 

Such a perturbation model is shown in 
sketch (v) . The sides of the pertur- 
bation body are flat rectangular sur- 
faces which are deflected. The wing 
end plates are not deflected in this 
case. A more efficient combination 
could probably be obtained by deflect- 
ing them the proper amount. 4 



■^he fins need not be on the wing. In this case the calculations are 
simplified in the tip region and the model is completed by the addition of 
control surfaces. In a design, the end plates may be replaced by external 
stores which are contoured to follow the curvature of the streamlines 
induced by the fuselage. The external stores would then receive a thrust 
from the wing and fuselage interference similar to Models 6 and 7 • 
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Lift.- The lift is divided into two parts: that due to the fuselage 

contour and that due to the wing deflection. The lift and drag generated 
by the wing deflection are given by equation (31). The indirect lift 
induced by the fuselage contour on the wing is calculated by the same 
method used for the other cases. The general equation for the basic coni- 
cal flow is obtained from solution 1 of table. A, 13a of reference 8. The 
stresmwise perturbation velocity for this conical flow is given by equa- 
tion (9b) . The pressure induced on the top of the wing by the fuselage 
is then 
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The indirect lift on the wing is found by integrating equation (4l) over 
the wing surface. It must be noted that the pressure waves generated by 
the fuselage are reflected by the fin at the wing tips, resulting in a 
doubling up of the lift in this region. After performing the integration 
with respect to x, the lift in the two regions is 
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where z is defined in equation (6c), 
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Drag .- In this case, the fins do not contribute' to the drag because 
they are undeflected and of zero thickness. The deflected plates on the 
body sides have a drag contribution from both their own deflection and 
the interference of the wing on them. The direct drag is found by inte- 
gration of the pressure as given by equation (9c) over the surfaces on 
the body. 


-'Dv- 


'v-v 


5 g c 
b + c 



3-36 


8 2 c 
b + c 


(^3a) 


The drag arising indirectly from the interference of the wing on the body 
is found by integrating equation (2l) over the area on the body influenced 
by the wing, as indicated by the crosshatched area in the side view in 
sketch (v) . 


P _ Son 
°h-v b + c 


(43b) 


Model 9 - Positive and Negative Wedge Superimposed 
on Fuselage; Subsonic Leading-Edge Wing 






Description .- The wings with 
sonic leading edges considered so 
far are not so efficient as some 
other plan forms. In an attempt to 
ascertain the effect of wing plan 
form on the efficiency of these 
special models, a more efficient wing 
is considered. The combination con- 
sists of a contoured fuselage and a 
wing with subsonic leading edges and 
sonic trailing edges. (See sketch (w).) 

The trailing edge Is taken to be sonic 
in order to simplify the geometry and 
the calculations. The taper ratio of 
the wing is l/4. The streamwise length 
of the rectangular sections on the fuselage are only long enough so that 
their influence does not reach the leading edge or extend past the trailing 
edge. This streamwise length is then equal to the root chord of the wing 
times the taper ratio. 



Both the wing and wedge are assumed to be deflected in the opt imum 
ratio which is defined by equation (29) . 

Lift.- The lift generated by deflection of the wing is read directly 
from figure 10 of reference 14. For the range of aspect ratios considered 
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the lift induced "by the body sections is calculated by equation (k2a), 
where b and the area of the wing are adjusted to this configuration. 
The loss in lift generated by the wedge in the region of the wing tip 
Mach cone is neglected. The final result would probably not be c hang ed 
by taking this into account. 

Drag .- Drag contributions arise from the wing deflection, wing 
leading-edge thrust, wing interference on thg fuselage, and from the 
deflected sides of the fuselage. The drag component due to the wing 
deflection is simply aC-j^ . The leading- edge thrust is calculated in 
the usual manner (see, e.gT, p. 219 of ref. 15) . The drag of the side 
plates is found from equation (43) with the denominator adjusted to the 
proper area. The expression for the pressure acting on the body due to 
wing interference is found from number 6 of table A, 13a of reference 8, 
as 


u = Re 


E(n/i - m2) 



m 


where E (vl - m 2 ) is the complete elliptic integral of the second kind 
(see ref. 16) . The real part of equation (44) iB found through the use 
of equation (4) and the pressure on the body is. 


( C p)y»0 * 


2a •V x 2 - zB cot g A 

E(n/i - cot 2 A) n/z 2 (1 - cot 2 A) + x 2 cot 2 A 


( 45 ) 


where the coordinate system and the area of integration are shown in 

sketch (x). This expression is integrated 
graphically with respect to z. 
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Model 10 - Half-Cone Mounted 
Under a Swept Wing 

Description .- So far the shapes 
used to generate a pressure field on 
the wing have been flat plates. Any 
shape whatever may be used to generate 
this interfering pressure field. One 
such example is one half of a body of 
revolution mounted on the bottom of a 
wing. The flow is then axially sym- 
metric and the model can be readily 
analyzed. A convenient choice which 
is relatively easy to calculate is the 
model wherein one half of a cone is 
placed under a swept wing as shown in 
sketch (y) . The wing leading and 
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trailing edges and the wing tip are swept along the Mach line. The wing 
may go to angle of attack as in previous cases. The cone is assumed to 
be slender and to rotate with the wing so that the boundary conditions 
on the wing and cone are not violated . 

Lift.- The lift on the wing arising from the deflection of the wing 
is found from equation (30) . Equation (31) cannot be used because it has 
the negative lift contribution of the third plan form shown in sketch (o) . 
The lift for this wing is, after integration of the pressure over the plan 
form shown in the bottom view in sketch (y) , 


■%- 


_ 8a P 

" 3* L 


2(b+c) ^2bc+c z -c z 


be + c 2 


(46) 


The lift generated by the cone interaction on the wing is found by 
integrating the pressure induced on the wing by the cone. The linearized 
pressures on the wing are given by (see, e.g., ch. 4, ref. 15 ) 


n _ ( 2u y£+w a 

°P VU “u 2 


°P 


= - 6 2 


S 2 (x 2 -P 2 r 2 ) _ 2 COBll -i _x 
r 2 (3rJ 


(47) 


where 6 is the tangent of the cone half- angle. Upon integration of 
equation (47) over the bottom of the wing, the lift induced by the cone 
is, for 0 = 1, 


26 2 c 2 


J c-h 2bc+c 2 1_ 


3 & cosh' 


-1 1 _ 3® 26f _ 


1-8 - 1 + 


!> 


28+1) 


bCs+ajcosh - 1 s±i + - » a m n 


zf2h+c\ ^ 2b+c 


2b 


6b 


^c+b^ f 2b+cT , _ 2 — cosh ~ x ^±5. 


3c 2 ~ c 2 


and the aspect ratio is 


A = 


2b +c 


(48) 

(49) 
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Drag . - The drag of- the wing is expressed by equation (31b). The drag 
of the cone alone is. 


C D, 


Jt5 4 c 2 


2(2bc + c 2 ) 


2 cosh" 1 ^ - 1 + 6 2 
o 


! ) 


(50) 


However, a drag also arises because the wing induces a pressure field 
around the cone.' In a first approximation, this drag may be taken as 
the frontal area of the half-cone times the pressure induced by the wing 
along its center line, that is, from equation (20) at y = z = 0. 



2x8 2 c 
2b + c 


(5D 


If a is taken as some product of 5, the optimum ratio of the wing 
angle of attack to the cone half-angle is found by equation (29) and the 
optimum drag factor from equation (28) . 


DISCUSSION OF RESULTS 


A sketch of each of the models analyzed, in the previous sections is 
shown in figure 1. The lift and drag coefficients and the parameter cr 
as predicted by linear theory for the models without fuselages are shown 
in figures 2 and 3. The values predicted for the models with fuselages 
are shown in figures ^ through 7« These results may be used by a designer 
interested in an aircraft utilizing interference effects between surfaces 
resembling those studied herein. 

The purpose in considering the models shown in figure 1 is to study 
the effect of the shape, location, and deflection of the various surfaces 
on the lifting efficiency of the airfoil system. 5 The question is whether 
indirect lift devices are worth considering. 

The lifting characteristics of the systems shown in figure 1 can best 
be compared by considering the drag factor Oj^/pCl 2 (see eq. (2)) as dis- 
cussed earlier. The drag factor for the models analyzed is shown in 
figures 8, 9 , 10, and 11 as a function of reduced aspect ratio J3A. The 
curve for the triangular wing is shown in the figures as a reference curve. 
It Is well known that when the flow velocity perpendicular to the leading 
edge is supersonic, (3A > 4, the lift-curve slope is 4/j3 and the drag 
factor is O.25. Also shown as -a reference curve is the theoretical mini- 
mum value of the drag factor for slender planar wings. This minimum or 
optimum curve serves to indicate whether the airfoil system being con- 
sidered is more efficient than is theoretically possible with a planar 
system. 

5 The effect on performance when the airplane is flying at an off- 
design Mach number is not considered. ~~ 
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It is Immediately obvious from figures 8 (a) and 9(a) that, when the 
drag factor is based on the plan form plus fin area ; the models are not 
nearly so efficient as a good planar wing. The figures 8 (a) and 9(a) 
apply to the physical cases therein the fins are added to the aircraft. 

The more promising situation is when the coefficients are based on the 
plan-form area (figs. 8 (b) and 9 (t>))> that is, when the vertical surfaces 
are assumed to be usable for other purposes , such as lateral control. 

The remainder of the discussion will be restricted to the cases when the 
coefficients are calculated using the plan-form area as the reference 
area. In other words, it is assumed that the vertical surfaces are a 
necessary part of the aircraft for reasons other than lift augmentation. 

The curves of the drag factor (figs. 8 (b) and 9(t>)) for the models 
when all the surfaces are deflected in the optimum ratio (solid curves) 
are noted to be the lowest of the results calculated. At low values of 
reduced aspect ratio these models are superior to the optimum planar wing 
expressed by the theoretical minimum curve. These results point out two 
important conclusions: (l) The lifting efficiency of aircraft can in 

certain cases be improved over planar systems by indirect lift devices . 

(2) The largest improvement is achieved if all of the available surfaces 
of the airfoil system are deflected in the correct ratio (see, e.g., 
eqs . ( 28 ) and ( 29 )). 

Various modificiations were made to the simple airfoil system. 

Model 1, and each represents an improvement in the lifting efficiency or 
a reduction in the drag factor as shown in figure 10. The addition of 
fins (undeflected) to the reverse swept wing results in a notable improve- 
ment in the lifting efficiency. A further gain is achieved by deflecti ng 
the fins and wing properly, and by increasing the lift by removing that 
portion of the wing which is carrying a negative indirect lift force from 
the fin, as discussed in the analysis of Models 2 and 7* This illustrates 
the importance of mailing certain that the entire surface of the airfoil 
system works positively. Another reduction in the drag factor tells one 
to contour the fuselage, even at the expense of losing the overlap lift 
in the center of the wing. Models 6 and 7* All of the surfaces of Model 7 
have been cut and deflected to produce the best configuration shown in 
this figure. 

The essential requirements of the wing of the airfoil systems studied 
are that it lie in the region of influence of the fin and that it serve as 
a barrier between the opposite pressure fields induced by the upper and 
lower fins. These considerations and the fact that the airfoil systems 
must be reasonably simple to analyze restrict the wing plan form somewhat. 
Therefore, with three exceptions, the wings studied in this report are of 
constant chord with leading and trailing edges swept along the Mach lines. 
It is to be noted in figure 8 , that the drag factor of this sonic leading- 
edge wing is not near the theoretical minimum for planar wings. The prob- 
lem treated is one of taking a wing which is not too efficient and trying 
to make it better than something which is already good. An alternative 
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problem Is to take a good wing and try to improve it by addition of body 
contouring and/or fins. This is attempted in. the design of Model 9, but 
the choice of the wing and body contour is not good. In order to facili- 
tate the analysis of the model the body contour had to be taken small 
enough so that the pressures induced by it did not bleed over the leading 
edge. Since the taper ratio of this wing is 1/4, only 2/5 of the wing 
surface felt the influence of the body contour. As a result, the wing 
alone and the combination are about equally efficient. Both results fall 
on the single curve shown in figure 9* The conclusion drawn from this is 
that when a device which is not too efficient is added to one which is 
efficient, the optimum for the combination is only slightly better than 
the efficient part by itself. The largest improvement is realized when 
the fin and wing are about equally efficient. A combination using a 
good wing approaching the theoretical minimum which could also utilize 
indirect lift to a great extent is the desired combination. 



One such airfoil system which can 
be analyzed theoretically is shown in 
sketch (z). The modjel consists of an 
arrowhead wing with subsonic leading edges 
and sonic trailing edges with the indirect 
lift induced by flat vertical sections on 
the fuselage sides which are juBt high 
enough so that their upper, and lower ends 
do not influence the wing. The extreme 
height of these flat sections enables one 
to find the pressures on the wing quite 
simply from linearized conical flow theory® 


6 The conical flow solution (provided by Doris Cohen of theJhnes Aeronautical Laboratory staff) 
for the axial velocity Induced on the wing by the fuselage wedge is 


g - Ref tan ' 1 + * 2 jjjL 

U « L 'Z l-m 2 4 m 2 -(y/x) 2 Ej 


where K and E are complete elliptic integrals of the first and second hind vith modulus s/l-m 2 . 
The lift and drag are then found as 


sfiSr * am 


[m 3 + ~= sin * 1 a + — gg_ ] "l 
L s/l-m 2 J 


CDv * v tm: 


88 c (*y l-m 2 - 1+m) 26 2 /l-m f m 2 ]! 2 * '/l-m 2 , 

I + -ar^TwlT - + 2 * E J 


texm /k 
1 +m \E 
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The egressions for ^ and ^ become more complicated for reduced aspect ratios below 

6 (m = 0.6) . 
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but makes the proportions of the model 
impractical. The drag factor for the 
model when the wing is undeflected 
(a = 0 ) and the wedge is present (8 / 0 ) 
and when the wing is deflected (a ^ 0 ) 
in combination with the wedge (8 / 0 ) is 
shown in sketch (aa) . The results for 
the arrowhead wing can be found from ref- 
erence 8 . A sizable improvement over the 
arrowhead wing is achieved by deflecting 
the fuselage sides in combination with 
the wing. However, thus far, the analy- 
sis can be carried out only for rather 
high fuselages and a model of comparable 
efficiency which is of more practical 
proportions is yet to be desired . 

The several different fin configurations which were studied permit a 
limited discussion of the effect of fin shape on the lifting efficiency. 

It is seen from figure 8 (b) that the airfoil combinations with unswept 
leading-edge fins. Models 1 and 2, are more efficient than when the leading 
edge is swept along the Mach line as with Model 3 • This is true for all 
values of reduced aspect ratios for the curved-trailing-edge wing plan 
form. Model 2, but only for the lower values for the constant -chord wing, 
number 1. It is interesting to note that although the drag of the fin 
shapes for Models 1 and 3 are the same, the pressure shadow on the wing 
is very different. 

For wings of reduced aspect ratio above 10, it may be advantageous 
to go to taller fin shapes such as those of Model 4. The rectangular 
section of the fin increases the region on the wing where the pressure 
distribution is of the two-dimensional type. The extent of this region 
increases roughly as the square of the fin height (see eq. ( 23 )), resulting 
in an efficient model at large aspect ratios . The drag factor for Model 4 
is not presented since the values are off scale for values of reduced 
aspect ratio below about 8 . 

If the fins were to lean outward a .slight amount, they would not only 
generate indirect lift on the wing, but carry a component themselves. A 
simple problem to investigate is the limiting case shown in sketch (m) as 
Model 5 . This case is the same as twisting the wing and does not permit 
use of vertical surfaces already on the airplane. The results of Model 5 
may be applied to wing-tip control surfaces . 

The vertical surfaces on the aircraft might be modified to receive 
a push from the pressure field induced by the wing or another vertical 
surface. This is illustrated by the wing-body combinations 6 and 7 
where the fuselage is contoured to receive a thrust from the fin and wing 
pressure fields. Typical values of the interference of the fin on the 
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body are shown in figure 5 and the resulting thrust is shown in figure 6. 
This principle could be applied to other types of fins and wings which 
appear more structurally feasible . 

Model 8, sketch (v), is designed with the thought that it would be 
applicable to conventional aircraft as shown in sketch (q) . The fuselage 
is contoured to generate indirect lift on the wing and the fins serve as 
end plates. The control of the airplane is achieved by simply moving the 
fins in the proper combination to produce turning, rolling, or pitching 
motions. Roll and pitch control is also obtained by deflection of the 
fuselage sides. The parameters calculated for fin deflection only may 
be used in the design of this type control system. An advantage worth 
noting is that the lift can be changed without first pitching the airplane 
or missile. 

The drag factor for the indirect lift from the half-cone is shown 
for several cone angles in figure 11. (The lower solid curve is the drag 
factor for 5 =,15° and the wing deflected the correct amount.) It may 
be reasoned that the half-cone body is the entire fuselage rather than a 
perturbation body. The pressure drag of the half-cone would then appear 
in equation (2) as part of C-^ instead of part of Cp. . The dotted 

line in figure 11 is the drag factor for Model 10 when the half-cone body 
is taken as the entire fuselage. A comparison of the results in figure 11 
with those in figures 8 and 9 shows that the efficiency of the half-cone 
model is not so good as seme of the other airfoil systems analyzed. This 
can be attributed to the fact that the pressure shadow induced by a body 
of revolution Is not so intense as that from a flat surface; and, the wing 
used Is not so efficient as the one used on the other models. 

It might also be asked whether the gain from the end plates at super- 
sonic speeds Is greater than at subsonic speeds. The results in refer- 
ence 2 show that at subsonic speeds the lift is increased enough to offset 
the added frictional drag of the fins. However, it is also pointed out 
that approximately the same gain could be realized by laying the end 
plates horizontally as extensions to the wing span. At supersonic speeds 
the percentage gain realized when the end plates are used as extensions 
to the wing span is not so high as that possible when the same area Is 
used as end plates . 

The models discussed were analyzed by linear theory in which it is 
assumed that the direction of the Mach waves throughout the disturbed 
flow field is the same aB in the free stream. The plan-form geometry was 
designed around this assumption which, to a first order for thin wings 
at small angles of attack, is permissible. However, in these special 
designs the location of particular surfaces must be such that they can 
benefit from the influence of another surface. The boundaries of these 
Influence regions or pressure fields lie along shock or Mach waves from 
the deflected surfaces. Hence, if the local wave angles differ greatly 
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from the free- stream, angles and the shape or location of the surfaces is 
not adjusted accordingly, a part of the desired gain may he lost. S imi lar 
considerations apply to flight at an off -design Mach number. 

Two examples are shown in figures 12 and 13 so that the comparison 
can be made between the wave pattern calculated by linear theory and the 
more exact wave patterns computed by the shock- expansion method. The 
model shown in figure 12 is similar to Model 4. It is to be noted that 
the compression regions lie ahead of and the expansion fan behind the 
linearized locations. Also note that the upper and lower edges in the 
more exact theory do not match. 

A bottom view of the linearized and the more exact wave pattern for 
a 10 ° half- cone mounted on the bottom of a swept wing is shown in fig- 
ure 13 . The more exact wave pattern was found by the method of character- 
istics. The line in the pressure field along which the lifting pressures 
go to zero is indicated as the Cp = 0 line. It is to be noted that the 
linear wing plan form lies within the lifting pressure field. 


CONCLUDING EEMARKS 


It is hoped that the foregoing analysis will provide some interesting 
and instructive examples of the type of interference effects that may be 
utilized at supersonic speeds. While no far-reaching conclusions should 
be drawn from the examples given, it is clear that in some cases the lift- 
drag ratio can be improved noticeably if interference effects between 
existing vertical surfaces and the wing are UBed advantageously. Several 
of the models studied are more efficient at low aspect ratios than is 
theoretically possible with planar wings. A configuration of comparable 
efficiency with a very high narrow fuselage is found for the higher values 
of reduced aspect ratio. The efficiency of the whole airplane can, in 
certain cases, be improved by indirect lift devices which can also be used 
for control purposes. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 8 , 1955 

REFERENCES 


1. Reid, Elliott G.: The Effects of Shielding the Tips of Airfoils. 

NACA Rep. 201, 1925 . 

2. Hemke, Paul E.: Drag of Wings with End Plates. NACA Rep. 267 > 1927- 



_ n n^TP-T T> ^ TTmTErUTUV T 


NACA RM A55L08 


40 


3. Scholkemeier, F. W.: Flugel mit Endscheiben. Ministry of Supply, 

Volkenrode Monograph, Reps, and Trans. 54l , 1948. (Also published 
as Ministry of Supply, Volkenrode Monograph GDC/2^72) 

Gorgui, M. A.: The Effect of De3ta Vanes on Supersonic Wings. The 

Aeronautical Quarterly, vol. V, pt. 4, Nov. 1954, pp. 251-279. 

5. Clements, Harry R.: Canted Adjustable Sid Plates for the Control of 

Drag. Aeronautical Engineering Review, vol. 14, no. 7, July 1955, 
pp. 40-44. 

6. Ferri, Antonio, Clarke, Joseph H., and Casaccio, Anthony: Drag 

Reduction in Lifting Systems by Advantageous Use of Interference. 
PIBAL Rep. 272, Polytechnic Institute of Brooklyn, Dept, of Aero- 
nautical Engineering and Applied Mechanics, May 1955* 

7. Lagerstrom, P. A.: Linearized Supersonic Theory of Conical Wings. 

NACA TN 1685, 1950. 

8. Jones, R. T., and Cohen, Doris: Aerodynamics of Wings at High Speeds. 

(Sec. A of Aerodynamic Components of Aircraft at High Speeds. 

vol. VII of High Speed Aerodynamics and Jet Propulsion, A. F. Donovan 

and H. R. Lawrence, ed., Princeton Univ. Press, In Press.) 

9. Cohen, Doris, and Friedman, Morris D.: Theoretical Investigation of 

the Supersonic Lift and Drag of Thin, Sweptback Wings With Increased 
Sweep Near the Root. NACA TN 2959, 1953. 

10. Harmon, Sidney M., and, Jeffreys, Isabella: Theoretical Lift and 

Damping in Roll of Thin Wings with Arbitrary Sweep and Taper at 
Supersonic Speeds. - Supersonic Leading and Trailing Edges. NACA 
TN 2114, 1950. 

11. Cohen, Doris: Formulas for the Supersonic Loading, Lift, and Drag 

of Flat Swept-back Wings with Leading Edges Behind the Mach Lines. 
NACA Rep. 1050, 1951. 

12. Hayes, Wallace D.: Reversed Flow Theorems in Supersonic Aerodynamics. 

North American Aviation Inc. Rep. No.“AL-755> Aug. 20, 1948. 

13. Heaslet, Max. A., and Spreiter, John R.f Reciprocity Relations in 

Aerodynamics. NACA Rep. 1119 , 1953* 

14. Malvestuto, Frank S., Jr., Margolis, Kenneth, and RIbner, Herbert S.: 

Theoretical Lift and Damping in Roll at Supersonic Speeds of Thin 
Sweptback Tapered Wings with Streamwise Tips, Subsonic Leading 
Edges, and Supersonic Trailing Edges. NACA Rep. 97^, 1950 • 





NACA RM A55L08 


41 


15. Heaslet , Max. A., and Lomax, Harvard: Supersonic and Transonic Small 

Perturbation Theory. Sec. D of General Theory of High Speed Aero- 
dynamics, W. R. Sears, ed., pp. 122-344. Vol. VI of High Speed 
Aerodynamics and Jet Propulsion. Princeton Univ. Press, 1954. 

16. Byrd, Paul F., and Friedman, Morris D.: Handbook of Elliptic 

Integrals for Engineers and Physicists. Vol. LXVTI of Grundlehren 
der Mathematischen Wissenschaften in Einzeldarstellungen. Berlin, 
Springer-Verlag, 1954. 





42 


CQHEIBEMW-fttnwfc 


NAjCA EM A55L08 


.. r niii¥iiTiwrnifiTii 



3 














Only vertical surfaces 

UU IIUVIUU , 

All surfaces deflected 

in optimum ratio . 



M'ACA BM A55LO8 








iWfi 






Gniy Vertical surfaces 

deflected . 

AH surfaces deflected 

in optimum ratio . 





NACA RM A55L08 




















uuip 



NACA JIM A55L08 




(a) Coefficients "based on plan-form plus fin area. 


Figure 8.- Drag factor for several airfoil combinations without fuselages. 
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Figure 13.- Plan view of linearized and true wave patterns for a half cone mounted under a Bwept 

wing; Moo » 2.0. 




